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Abstract 

We study the formation of knots on a macroscopic ball-chain, which is shaken on a horizontal 
plate at 12 times the acceleration of gravity. We find that above a certain critical length, the 
knotting probability is independent of chain length, while the time to shake out a knot increases 
rapidly with chain length. The probability if finding a knot after a certain time is the result of the 
balance of these two processes. In particular, the knotting probability tends to a constant for long 
chains. 

PACS numbers: 



1 



I. INTRODUCTION 



Knots are prevalent on most cables, chains, and strings being used in every day life or 
technology Most remarkably, knots appear to form spontaneously, as soon as strings are 
shaken, transported, or handled in any way, and thus are an unavoidable byproduct of their 
use. For example, [l] describe the dynamics of a ball chain that is suspended from an 
oscillating support. As soon as the chain dynamics become chaotic, the chain forms knots 
of various kinds. Yet in spite of a considerable amount of work on the importance of knots 
on the m oleeufa r sea, flfi, we a r e act aware of any S _,o study into the on gln of the 
prevalence of knots on macroscopic chains. 

In this paper we present model experiments on shaken ball chains that quantify the 
tendency for knot formation as function of chain length. By considering both knotting and 
unknotting events, we present a simple theory that explains the probability for the formation 
of knots after the chain has been shaken for a given amount of time. After the experiment, 
the ends of the chain lie flat on the plate, and there is a unique way to join the ends to form 
a closed curve. If this curve is topologically equivalent to a closed loop or "unknot" '\ there 
is no knot, otherwise we call the chain knotted. We made no distinction between different 
kinds of knots, however the simple trefoil knot [4] was by far the most common. 

Of course, it is precisely the topological stability of knots that lies at the root of the 
phenomenon: once a knot is created, it cannot disappear, except when it falls out at the end 
of the chain. Using a setup very similar to ours, the pioneering study [5j] investigated the 
lifetime of a simple trefoil knot that was placed in the middle of the chain at the beginning 
of the experiment. The mean lifetime r of the knot was found to increase rapidly with chain 
length: 

t = t 3 (N-N ) 2 /D, (1) 

where N is the number of balls on a chain, N w 15 the size of a knot, and D a hopping 
rate. The knot was modeled by the three points of intersection of the chain, which perform 
random walks. From this assumption the constant r 3 = 0.056213, as well as the entire 
distribution of lifetimes was calculated. The experimentally determined hopping rate of 
D — 11 ± Is -1 was found to be close to the driving frequency of / = 13Hz; however, it has 
not been investigated systematically what sets this time scale. 

In Fig. we present a sketch of our experimental setup. A stiff solid plate, 50 cm 
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FIG. 1: Sketch of the experimental setup. 

in diameter, is vibrated vertically by an electrodynamic shaker. The plate's weight of 4 
kg is sufficiently large for the chains to have little effect on the motion of the plate. The 
plate is attached to a metal block mounted on a spider spring, to support the weight of 
the plate and to reduce lateral motion. The shaker is controlled by a computer, with an 
accelerometer providing feedback to make sure the plate motion is close to sinusoidal. The 
driving frequency / = 21Hz is dictated by the requirement of operating close to resonance, 
the dimensionless acceleration is T = Au 2 / g = 12 ± 0.2 as recorded by the accelerometer. 
Here A is the amplitude, u = 2irf the angular frequency, and g the acceleration of gravity. 
The top of the plate was machined to have a very shallow parabolic profile, of 5mm depth 
at the center of the plate. This amount of confinement was enough to always keep the 
chains near the center, without them ever feeling the edge of the plate. A digital camera 
was mounted above the plate, capable of taking up to 20 frames per second. 

Ball chains are an excellent model system to study knot formation |5j, in that they have 
little stiffness that would resist the formation of loops, yet considerable friction between the 
beads keeps knots from opening too easily. All our chains were cut from a single sample 
that had a bead diameter of 2Ri iea( i = 2.1 ± 0.05mm, with connecting rods that allowed 

for a maximum inter-bead spacing of 2.1 ± 0.05mm = 2i4ead, somewhat greater than that 

n 

of the earlier study (0.8Rbead) 5]. It seems reasonable to take the number of beads as the 
fundamental unit of length, since they offer the greatest resistance as two parts of the chain 
slide across each other. 
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FIG. 2: Probability for knotting after 30 seconds. Error bars are calculated from the variance of a 
binomial distribution. The full line is based on equation (J3J) below. 

II. EXPERIMENTAL RESULTS 

Our main aim is to understand and to quantify the tendency for knots to form sponta- 
neously once excited by shaking. To that end we shake chains of lengths between N=10 
and 500 by dropping them onto the vibrating plate, and inspecting them for knots after 30 
seconds. Each experiment was repeated 80 times, and the resulting probability of knotting 
calculated, as shown in Fig. |21 No knotting was ever observed for chains shorter than 
Nmin = 38. The probability then rises sharply to reach a plateau value of about P = 0.26. 

We expect that Fig. |2] can be understood from the interplay of knotting and unknotting 
events. To investigate this further, we followed the evolution of chains whose length lay in 
the transition region (N=73 and N=144) continuously for two hours. By taking 10 frames 
at 20 fps we were able to decide unambiguously whether a knot was present. This process 
was repeated every 5.5 seconds, which permitted us to compress and store the video images. 
The entire sequence was then examined manually for the number of knots, separated by a 
length of chain. As shown in Fig. El the presence of two knots is still quite unlikely for the 
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FIG. 3: The number of knots on chains of length N=73 and N=144 as a function of time. 



chain lengths shown here. In calculating the statistics of knotting and unknotting, knots 
were treated as independent. 

Some statistics as extracted from FigJ3]are collected in Table |U As expected, the lifetime 
of a knot increases considerably with length. Comparing to the results of [5], there is 
remarkably good agreement, using with the same value of the hopping rate D as in the 
original paper. Th,s is perhaps 8 u rp ris,„ g , considering the differences in chain propert.es, 
driving conditions and, in particular, in the way knots are introduced. In p, knots were 
introduced by hand in the center of the chain, while spontaneous knots tend to originate 
from the end of the chain, see Fig. 0] below. Note that our driving frequency was also 
somewhat larger than that of Q]. 

However the most remarkable observation is that the mean knotting time for the two 
chains is the same, although the chain lengths are quite different. This observation makes 
sense, since knots are produced by the ends of the chain, which have sufficient freedom of 
motion to wrap around the rest of the chain, as illustrated by the three knotting events shown 
in Fig. 0] Of course, there is certain minimum length that is required, which according to 
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Chain length [N] 



73 



144 



Chain length [cm] 30 60 

mean knot lifetime [s] 13 ± 5 102 ± 29 

r [s] [5] 17 85 

mean knotting time [s] 219 ± 12 219 ± 16 

< R g > [cm] (no knot) 7.18 ± 0.03 11.58 ± 0.06 

< R g > [cm] (knot) 6.39 ± 0.12 11.00 ± 0.09 



TABLE I: Some statistics of knotting and unknotting for two chains, N=73 and N=144. The last 
two lines report the mean radius of gyration. 




FIG. 4: Examples of knotting events, taking place at the end of the chain. 

Fig. El is N m i n = 38, about twice the minimum size Nq of a knot. Our data for N = 73 
indicates that the knotting rate rises very quickly to a plateau value in the chain length. We 
will assume that the knotting rate is actually constant once N is larger than N min . In Fig. 
EJ we plotted the cumulative distribution of knotting times for the shorter chain, as obtained 
from the data of Fig. El The result is very well fitted by an exponential distribution, using the 
mean knotting time from Table H] This is the distribution we will assume below for both 
knotting and unknotting events, although the unknotting distribution is more complicated 



As a final result, we report the mean radius of gyration < R g > obtained by taking the 
spatial average of all beads of a chain, which was done by image analysis using MATLAB. 
We then performed a temporal average, which we did separately for the knotted and the 
unknotted chain. As expected, < R g > is slightly larger for the unknotted state, but 
temporal fluctuations are prohibitively large to make this a useful indicator of knottedness, 
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FIG. 5: Cumulative number K of knotting events occurring after a shaking time of t or smaller. 

as we had hoped originally. The radius of gyration of course increases with chain length, 
but at a slower rate than the linear extension of the chain. 



III. THEORY AND DISCUSSION 

Armed with the above observations, we can attempt a simple theory for the probability 
of knots. We adopt a two-state description, in which the chain is either in an unknotted 
or in a knotted state, the probability of the latter being P. This can easily be generalized 
to allow for an arbitrary number of knots. Events are completely uncorrelated in time, 
so the distribution of (say) unknotting times is assumed exponential, while reality is more 
complicated |5j. Now if the average unknotting time is r as before, and the knotting time 
Tfc, we obtain the following rate equation for P: 

P=(l-P)/r k + P/r. (2) 

The first term is the rate of knotting events, which only take place if there is no knot. 
Conversely, the second term describes the rate of unknotting. Equation (J2J) is solved very 
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simply with initial condition P(0) = 0, giving 



P(t) 



1 



(3) 



1 + T- fe /r 



for the probability after a time t. Now all that is needed are the values of r and as 
function of chain length. 

For r we take ((TJ as found in [^j], which according to Table H] is consistent with our data. 
For the knotting time we take 



As discussed above, this is based on the idea that knots are formed at the end of the chain. 
Once a length N min , which gives sufficient freedom for knots to form, is exceeded, the rest 
of the chain no longer matters. By the same argument we also believe that the mean 
configuration of the chain, as measured qualitatively by the radius of gyration, is not of 
fundamental importance to calculate knotting probabilities. The constant rjf at ' has to be 
determined empirically. 

For long chain lengths, r becomes large and ® reduces to P = 1 — exp(— t/r^ sat ^), 
independent of chain length, as expected. From the asymptotic value of P = 0.26, taken 
from Fig. El we deduce = 100s. For simplicity, we also assume that N min is essentially the 
same quantity as No as identified by j^, the length for which a knot falls out immediately. 
The result, (J3J), is plotted as the solid line in Fig. |21 using (HJ,(HJ). The agreement is quite 
good, considering there is only one adjustable parameter. 

Ideally, even this adjustment could have been avoided, as could be taken directly from 
Table |U However, this gives about twice the value of T&, which would lead to a significant 
disagreement in the asymptote of P. One simplifying assumptions of our model was to 
disregard the distribution of knotting times. However we suspect that the main reason for 
this disagreement lies in the difficulty of preparing the chain in an unbiased fashion. At 
the beginning the ends of the chain are excited to a higher degree, leading to a significant 
increase of the knotting probability. 

There is a significant body of work that remains to be done. Firstly, one would like to 
check our theory in greater detail, by taking long time traces for a greater variety of chain 
lengths, and measuring knotting probabilities for periods other than 30 s. Secondly, one 




N <= N, 



N> N, 



nun 



(4) 
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would like to obtain a better description of the dynamics of the chain, and how it leads 
to knotting events. In particular, what sets the time scale for r and ? This could be 
addressed in part by changing the driving characteristics of the plate, as well as the amount 
of confinement of the chain. Preliminary experiments, performed in a container with side 
walls, have in fact shown that boundary effects play a significant role for the knotting 
probability. Finally, we have not considered the probability for different types of knots, 
which are observed if the knotting time increases. 
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